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ABSTRACT 

We report an updated analysis of the gamma-ray source AGL J2241+4454 that was 
detected as a brief two-day flare in 2010 by the AGILE satellite. The high-energy 
emission of AGL J2241+4454 has been attributed to the binary system HD 215227, 
which consists of a Be star being orbited by a black hole making it the first known Be- 
black hole binary system. We have analyzed the AGILE data and find a gamma-ray 
flux of (1.8±0.7) x 1CU 6 ph crrU 2 s” 1 , in agreement with the initial report. Additionally, 
we examined data from the Fermi LAT over several time intervals including the two 
day flare, the folded orbital phase, and the entire mission (~6-years). We do not detect 
AGL J2241+4454 over any of these time periods with Fermi and find upper limits 
of 1.1 x 10 -7 ph cm -2 s _1 and 5.2 x 10 -10 ph cm -2 s _1 for the flare and the full 
mission, respectively. We conclude that the HD 215227 Be-black hole binary is not a 
true gamma-ray binary as previous speculated. While analyzing the Fermi data of the 
AGL J2241+4454 region, we discovered a previously unknown gamma-ray source with 
average flux of (13.56 ± 0.02) x 1CD 8 ph cm~ 2 s _1 that is highly variable on monthly 
timescales. We associate this emission with the known quasar 87GB 215950.2+503417. 

Key words: binaries: general - gamma-rays: stars - stars: black holes - stars: 
emission-line, Be - stars: black holes - quasars: general 


1 INTRODUCTION 

Microquasars are typically binary systems consisting of a 
normal (non-degenerate) star from which material is pulled 
into an accretion disk around a black hole (BH) compan¬ 
ion, though in some c ases the accr eting object may be a 
neutron star (NS; see Izhan 3 [ml for a review of micro¬ 
quasars). These objects may be considered as scaled down 
versions of extra-galactic quasars as they appear to share 
similar properties. They emit at a wide range of frequen¬ 
cies from radio through gamma-ray, where the emission at 
different energ ies is often correlated thro ugh disk-jet in¬ 
teractions (iRemillard fe McClintockl 120061 ). The accretion 
disk in microquasars c an be fed from Roche-lo be o ver- 
flow of the companion dMcClintock fc Remillardl 1 19861 ) or 
from dense the stellar w inds of a massive secondary star 
(Ivan Kerkwiik et al.lll993) . 

The objects Cyg X-l and Cyg X-3 are known micro¬ 
quasars and high-mass X-ray binaries (HMXB). While Cyg 
X-l is known to harbor a 15 Mq BH (lOrosz et al 1l201l|) , the 
nature of the compact object in Cyg X-3 is remains open, 
though some evidence suggests a 2 - 5 M a BH companion 
dShrader et all 120101 ; IZdziarski et all 12013 ). Both systems 
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have been shown to emit gamma-r ay flares dBulearclli et aid 
l20ld . l201ll ; iBodaehee et all 120131), with some evidence for 
orbital modulation in Cyg X-3 dAbdo et all l2009bl) . The 
gamma-ray emission is associated with the soft X-ray state 
in Cyg X -3 dBulearelli et al. 20121 ), the hard X-ray state in 
Cyg X-l dMalvshev et al.ll2013l ). with heightened radio emis¬ 
sion in both. 


Another class of objects are the gamma-ray binaries of 
which five are known: LS I +61°303, LS 5039, PSR B1259- 
63, HESS J0632+057, and 1FGL J1018.6-5856 dDubusI 
120131) . These objects exhibit high-energy spectra that peak 
at GeV and TeV energies, in contrast to X-ray binaries 
whose spectra at keV energies. It is suspected that the 
compact companions of gamma-ray binary systems are NS 
and that the gamma-ray emission arise s from the colliding 
pulsar and stellar winds (|Pubusll2013l ). However, gamma- 
ray emision may also arise from the interaction between a 
clumpy stellar wind and the jet produced by a microquasar 
dOwocki et al.l[2009l ). 

In 2010, the AGILE satellite dTavani et aid 120091) de¬ 
tected a gamma-ray flare with a flux of 1.5 x 10~ 6 
ph cm" 2 s _1 (>100 MeV) with a significance greater than 
5 a from a prev i ously unknown source, AGL J2241+4454 
dLucarelli et aid 120101 ). The location of the flare co¬ 
incides with a known Be star, HP 215227 (=MWC 


© 2014 RAS 












































2 Michael J. Alexander and M. Virginia McSwain 


656; iMerrill fe Burweli 11943h . whose binary nature was 
elicited only after the discovery of the gamma-ray 
flare ([Williams et all I2010h . This began speculation that 
AGL ,12241+4454 could be the sixth identihed gamma-ray 
binary. I Mori et al.l d2013l l performed an analysis of Fermi 
data and found upper limits (90% C.L.) of 9.4 x 10” 10 
ph cm -2 s” 1 for the 3.5-year average flux and 7.2 x 10” 8 
ph cm -2 s” 1 for the flare. They also investigated orbital 
modulation by folding the data over the period of the bi¬ 
nary and found upper limits for each phase bin in the 
neighborhood of 10” 8 ph cm -2 s” 1 . So far, observations 
from Fermi have been unable to confirm the detection of 
AGL J2241+4454. 

Recentlv. lCasares et al.l (1201414 found the Be star’s com¬ 
panion to be a stellar-mass BH, making it the first Be 
star in a black hole binary (BHB) system and b ringin g 
the total number of confirmed BHBs to 25 dZhan 5 120131) . 
The secondary Be star has been classified as Bl.5-2 III 
dCasares et al.ll2014l ') , whil e a previous analysis fo und a spec¬ 
tral type of B3 IVne+sh dWilliams et aid 201Cf) . The mea¬ 
sured orbital period of 60.37 days and Be star mass of 10 -16 
Mg yield a BH mass of 3.8-6.9 Mg dCasares et al.ll2014lh 
HD 215227 has also been identifi e d as a source of X-ray 
emission by iMunar-Adrover et al.l d2014lh , who find a non- 
thermal X-ray luminosity Lx = 3.7 x 10 31 erg s” 1 corre¬ 
sponding to 3.1 x 10” 8 Lb dd for a 3.8-6.9 Mq BH, consistent 
with t he previously f ound upper limit of Lx < 1.6 x 10 7 Le<m 
dCasares et al.l 120141 ). This low luminosity is in stark con¬ 
trast to Cyg X -l, which has an Lx ~ 10 37 erg s” 1 and 
LEdd ~ 10” 2 dGallo et aid l20l5 l. This suggests that the 
system was in a quiescent state at the time of the X-ray 
observations. While there is a strong contrast at X-ray ener¬ 
gies, neither Cyg X-l nor AGL J2241+4454 show evidence 
for sustained gamma-ray emission, bu t have both be en de¬ 
tected as gamma-ray flares by AGILE dLucarelli et al.ll~201Cl : 
iBulgarelli et alJl20inh . Three separate AGILE detections of 
Cyg X-l we re confirmed as low-si gnificance counterparts in 
Fermi data l|Bodaghee et al.ll20 13h The highest significance 
detections were more likely to be found during transitional 
stages, perhaps coinciding with jet formation or destruction. 


The low X-ray and radio flux from HD 215227 
places near the low-mass X-ray binary (LMXB) A 0620- 
003 (Lx =~ 10 30 ' 5 er g s -i. | Garcia et al.l 120011 1 

on the well-formed radio /X-ray luminosity relation 
dMunar-Adrover et alJ 2014 1. T he X-ray luminos ity o f 
HD 215227 (~ W~ s LEdd', IMunar-Adrover et al.l l2014l l 
is also found to be similar to the median value of sam- 
ple of eight quiescen t systems (Lx = 5.5 x 10 ~ 7 LEdd; 

I Reynolds et al.l l2014l l. In this quiescent state, the X-ray 
luminosity of these may be dominated by emission from an 
advection-domi nated accretion flow (ADAF) rather than 
jet interactions ([Reynolds et al.ll2014l f. 

In light of new evidence for AGL J2241+4454 as both 
a HMXB and a BHB, we present an updated analysis of 
gamma-ray emission towards AGL J2241+4454 including 6- 
years worth of Fermi data. We aim to determine whether or 
not AGL J2241+4454 belongs in the rare class of gamma- 
ray binaries. Section [2] presents the details of our analysis 
of AGL J2241+4454 using AGILE data, Section [3] includes 
updated analysis of ~6 years of Fermi data, Section 0 ] high¬ 
lights the detection of a newly identified gamma-ray point 


source in the nearby field, and Section [5] discusses and sum¬ 
marizes our results. 


2 AGILE ANALYSIS OF AGL J2241+4454 

The A GILE satellite mission launched in 2007 by the Italian 
Space Agenc y to survey the sky from hard X-rays through 
gamma-rays (iTavani et alj|2009f ). The main instrument for 
detecting gamma-rays is the Gamma-Ray Detecting Imager 
(GRID), which is sensitive in the energy range from 30 MeV 
- 50 GeV. The GRID PSF containment radius is 1-2° at 
300 MeV with an error box radius of 6-20' depending on 
source intensity and spectral shape. One of the main science 
drivers is rapid response to gamma-ray bursts and other 
high energy transient phenomena, so AGILE was designed 
for high precision timing and is capable of detecting strong 
(> 10” 6 ph cm -2 s _1 ) point sources with high significance 
in j ust 1 -2 days of o bserva tion. 

iLucarelli et al.l (j2010 i) reported the gamma-ray flare 
AGL J2241+4454 at the position (l = 100.0, b = -12.2) 
± 0?6 (stat.) ± 0?1 (syst.) lasting from 25 July 2010 01:00 
UTC through 26 July 2010 23:30 UTC. The measured flux 
during this time was 1.5 x 10” 6 ph cm -2 s” 1 at 5<r. To 
analyze this source we downloaded data covering that time 
period from the ASI Data Science Center websitcQ. Data 
reduction followed the GRID analysis manual for using the 
latest version of the AGILE software (v5.0). 

The position, flux, and spectral index of 
AGL J2241+4454 were set to the initial discovery values 
and allowed to vary during the maximum likelihood fitting. 
The resulting fit yield a position of (Z = 100?1,6 = 12?2), 
a flux of (1.2 ± 0.6) x 10” 6 ph cm -2 s _1 , and a spectral 
index of —1.3 ± 0.4 with a significance of 3<r. Because of 
the low overall photon count (7 ±4), we also performed the 
analysis with the spectral index fixed to a value of —2 to 
reduce the number of free parameters. This yields a flux of 
(1.8 ±0.7) x 10” 6 ph cm” 2 s _1 and a total photon count of 
11 ± 5 also at 3 <t, in agree ment with the initially reported 
value dLucarelli et all 1201 Oi l. 

We also investigated any spectral dependence by per¬ 
forming the same likelihood analysis on the different energy 
ranges available from AGILE. Figure |T] shows the measured 
fluxes broken into several energy bands: 50-100 MeV, 100- 
400 MeV, 400 MeV-1 GeV, 1 -3 GeV, 3-10 GeV, and 10-50 
GeV. No indivual band is detected above 3cr and all points 
are represented as upper limits. The solid black line shows 
a power-law slope of V = —1.3 while the dotted line has a 
slope of r = —2.0. These lines are representative and not 
fits as there are no solid detections in the individual energy 
bands. 


3 UPDATED Fermi ANALYSIS OF 
AGL J2241+4454 

The Fermi gamma-ray satellite was launched in 2008 by 
NASA to observe the high-energy sky. On-board are sev- 
eral instruments in cluding the Large Area Telescope (LAT; 
I Atwood et al.l[2009l) . The LAT has a 2.4 sr field of view and 

1 http://agile.asdc.asi.it/ 
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Figure 1. Gamma-ray flux upper limits measured by AGILE 
for different energy ranges. The horizontal grey lines indicate the 
bandwidth of each energy range. The solid black line shows a 
representative power-law with a slope of T = —1.3, while the 
dotted line is a power-law of slope T = —2.0. 


covers an energy range from 30 Mev to 300 GeV. We down¬ 
loaded and analyzed the Pass 7 Reprocessed data includ¬ 
ing photons within a 20° radius around AGL J2241+4454 
and energies from 100 MeV to 300 GeV. All analyses 
employed the Fermi SCIENCE TOOLS v9r32p5 package as 
well as the latest Galactic diffuse model (gll_iem_v05), 
isotropic spectral template, instrument response function, 
and point source catalog to which we added a source 
at the location of AGL J2241+4454. Followi n g previous 
studies of other ga mma-ray binaries l|Abdo et alJ [2009a; 
lHadasch et al.l [2012h . we modeled the source spectrum as 
a power law plus an exponential cutoff of the form, 



where No is the normalization, Eq is the scale factor, F is 
the power-law index, and E cuto ff is the high-energy cutoff. 

To analyze the flare, we restricted the Fermi good 
time intervals (GTIs) to the 46.5 hr duration of the flare 
llLucarelli et al.l I201CI ). Of this time, Fermi pointing and 
quality cuts limit the useful observations to 16.8 hrs of data. 
Our analysis shows no detection of the flare by Fermi , and 
we derive a 68% (la) confidence upper l imit for the flare of 
1.1 x 10 -7 ph cm -2 s _1 , consistent with lMori et alJ d2013h . 
Given the strength of the flare detected by AGILE (~ 10 -6 
ph cm -2 s _1 ) it is somewhat surprising that it was not de¬ 
tected, however, the low fraction of Fermi observations of 
the AGL J2241+4454 flare (~ 36%) may not have been 
enough time for the source to rise above the background 
flux. 

We performed a similar analysis for the ~6 years of 
Fermi LAT observations to look long-term gamma-ray emis¬ 
sion and find no evidence for sustained emission. We set a la 
upper limit to the average flux of 5.2 x 10 -10 ph cm -2 s^ 1 , 
nearly a factor of two lower than the previous upper limit, 
which likely owes to the additional data available. 


2 http://fermi.gsfc.nasa.gov/ssc/ 



Figure 2. Photon count map for the 6-year dataset (left panel) 
and the residual map (right panel). The 1° cyan circles show the 
location of AGL J2241+4454 while the dark blue circles mark the 
previously unknown source 3FGL J2201.7+5047. 
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Figure 3. Phase-binned gamma-ray fluxes for AGL J2241+4454; 
all points are la upper limits. The estimated AGILE flux for 
AGL J2241+4454 (1.5 X 10 -6 ph cm -2 s -1 ) is two orders of 
magnitude larger than the 6.5 yr average phase binned limits. 



Figure [2] shows the photon counts map on the left and 
the residual (counts minus model) map on the right. There is 
no obvious source of gamma-ray emission at the location of 
AGL J2241+4454 above the level of the bright diffuse emis¬ 
sion from the Galactic plane. There is an obvious residual 
leftover in the upper right of the image that is a previously 
unknown gamma-ray source (Section [3} • 

K nown g amm a-ray binaries such as LS I +61°303 
d Abdo et al.l l2009al ) show strong gamma-ray flux modula¬ 
tion over the orbit, so we searched for evidence of this from 
AGL J2241+4454. We divided the data into time bins of 
6.037 days, one-tenth of the orbital period, and folded it 
over the full time range covered by the data us ing a reference 
time of HJD = 2,453,243.7 as phase cj) = 0.0 dCasares et al.l 
120141) . Figure [3] shows the results of the phase-binned analy¬ 
sis. There are no signficant detections at any orbital phase, 
and all the points shown are upper limits. 

It is curious that AGL J2241+4454 was detected by 
AGILE and not by Fermi since the flux was measured 
above 10 -6 ph cm -2 s -1 . As mentioned earlier, one pos¬ 
sible explanation is that 16.8 hrs of on-source time for the 
LAT was insufficient despite the flare’s luminosity. To test 
this possibility, we ran a likelihood analysis on the source 
LS I +61 303, which has an average flux of 0.9 x 10 -6 
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ph cm -2 s _1 iHadasch et al.l[2012i) . comparable to that of 
AGL J2241+4454. LS I +61 303 was only observed for a to¬ 
tal of 4 hrs during the flare, so we chose a different time that 
matched the 16.8 hrs of AGL J2241+4454 observation. For 
this time period, we find a flux of 0.3 x 10~ 6 ph cm~ 2 s _1 
at the 4 a level. This situation is not entirely analogous to 
that of AGL J2241+4454 as LS I +61 303 emits consistently 
at ~ 10 -6 ph cm -2 s -1 , while AGL J2241+4454 may only 
have peaked at that flux for fraction of the flare.. However, it 
demonstrates that the LAT is able to detect bright sources 
in very short periods time. 

Another possible explanation is that the flare was not 
coincident with Fermi GTIs and the flux peak was not ac¬ 
tually observed. Figure [4] plots the AGILE and Fermi pho¬ 
ton arrival times versus the Fermi GTIs during the flare; 
the time begins at the reported flare start. Fermi photon 
counts (open diamonds) are summed for each individual 
Fermi GTI, which are marked by the grey shaded areas. 
The arithmetic median of the Fermi count rate is shown as 
the solid horizontal line, and simple la, 2a, and 3er standard 
deviations from the median are the dotted, dashed, and dot- 
dashed lines, respectively. AGILE photons are binned every 
60 seconds for better time resolution. The data represent all 
‘clean’ AGILE and Fermi source plus background photons 
within 5° of HD 215227. For the AGILE data, we followed 
the AGILE users’ guide to select the following additional cri¬ 
teria: event status G, Earth tolerance 5°, and orbital phase 
code 2. The Fermi data were extracted only after passing 
the photon files through the GTMKTIME step used for the 
primary data analysis. 

Figure 0] shows that the gamma-rays detected by AG¬ 
ILE arrived at the detector in the middle of the reported 
time limits with the peak flux near hour 22 (about 23:00 on 
25 July). However, of the 11 photons selected by the con¬ 
straints given above, just two landed within Fermi GTIs. 
The GTI immediately following the AGILE spike also shows 
a peak in the Fermi photon count, but did not reach the 
nominal 3<r limit. The figure cannot definitively demonstrate 
the presence/absence of a flare, however, it illustrates the 
possibility that the AGL J2241+4454 flare was not detected 
by Fermi because the telescope was simply not looking. If a 
flare did occur, then it must have lasted less than about six 
hours for the resulting flux not to have been seen in at least 
one additional GTI, including the maximum at 23.5 hours. 
A six-hour window translates to a light-travel distance of 
43 AU, which would be an upper limit on the size of the 
emitting region. 


4 THE NEW GAMMA-RAY SOURCE 
3FGL J2201.7+5047 

During our initial analysis of AGL J2241+4454, we found a 
bright object in the residual image (Figure [2j right panel). 
This source, which had no associated 2FGL source but ap¬ 
pears in the new catalog as 3FGL J2201.7+5047, we modeled 
as a log-parabola to the model file of the form, 

dN / 75 \ ~(a + Plog(E/E b )) 

m ‘ w ° (I) • < 2 > 

where Nq is the normalization, E b is the break energy, and 
a and /3 are power-law indices. We then followed the same 



Flare Elapsed Time (hours) 

Figure 4. Photon counts for Fermi (open diamonds) and AGILE 
(filled circles) during the flare. The Fermi data points represent 
non-background subtracted photon counts for the GTI (filled grey 
stripes), while the AGILE data was binned every 60 seconds. The 
horizontal lines mark the Fermi median (solid), and la (dotted), 
2 a (dashed), and 3a (dot-dashed) above the median. To is the 
reported flare start time of 01:00 UTC on 25 July 2010. 
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Figure 5. Fluxes (top) and significance values (bottom) for 
3FGL J2201.7+5047 with each time bin representing one month 
of Fermi observations. Diamonds represent valid (> 3cr) detec¬ 
tions while arrows indicate la upper limits. The vertical dotted 
and dashed lines indicate the time of the last data included for 
sources in the 2FGL and 3FGL source catalog, and the horizontal 
dot-dashed line in the bottom panel marks the 3a detection limit. 


binned likelihood analysis as before. The source was detected 
over the energy range 100 MeV - 300 GeV at 47cr with a 
flux of (13.56 + 0.02) x 10" 8 ph cm" 2 s _1 , a = 2.577+0.004, 
/3 = 0.098+0.002, and Eb = 505.2+0.7 GeV, while the Fermi 
catalog reports a 27er detection for 3FGL J2201.7+5047 of 
(0.219 ± 0.015) x 10~ 8 ph cm -2 s _1 , roughly two orders of 
magnitude lower. The main cause of the difference is most 
likely that the source remained in a high flux state (> 10 ' 
ph cm -2 s _1 ) for an 8-12 months after the data cutoff date 
for the 3FGL catalog (Figure [5]), which lead to a lower av¬ 
erage flux. Also, the 3FGL J2201.7+5047 was modeled as a 
power law using a smaller energy range (1-100 GeV) than 
for our analysis. These are the main factors that make the 
fluxes different, but not necessarily in conflict. Because many 
known gamma-ray sources are variable, we then divided the 
data into 30 day time bins to measure the changes in the 
high-energy flux. 
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Figure [5] shows the source flux since the start of Fermi 
science operations where each data point represents one 
month of observations. The dotted and dashed horizontal 
lines represent the time of last data included in the 2FGL 
(01 August 2010 at 01:17 UTC) and 3FGL (31 July 2012 
at 22:46 UTC) catalogs, respectively. During the time range 
included in the 2FGL catalog (dashed horizontal line), the 
detection significance never rose above 3<r in a given month, 
which is likely the reason the source was not identified. It 
appears to have brightened suddenly near the start of 2011 
and remained strong through the end of 2013. However,the 
last few months of data suggest that the object may be en¬ 
tering a low flux state. Continued monitoring is necessary 
to see if the gamma-ray flux continues to decline. 

We used the procedure GTFINDSRC to localize the posi¬ 
tion of the source to a = 330?412, 5 = 50?828 with an error 
circle radius of 0?016 (57”). The values agree well with the 
3FGL coordinates, a = 330? 434 ±0? 026, 5 = 50? 800 ±0? 025, 
where the errors are at the 68% confidence level. We re¬ 
port errors that are about two times smaller than the Fermi 
catalog, which may also be related to the additional high 
flux data in our analysis leading to a stronger detection 
and better localization. We searched the SIMBAD Astro¬ 
nomical Database^ around this position for possible pro¬ 
genitors and found a known blazar, 87GB 215950.2+503417 
(IGregorv fe Condon]Il99lh . 62" away. This association was 
also recently identified by the Fermi team in the 3FGL 
catalog. The quasar is also coincident with X-ray source 
1RXS 220140.8+504941 from the ROSAT All Sky catalog 
(IVoges et al.lll999l V). which has a count rate of (2.84+0.83) x 
10 -2 s -1 . Prior to its launch, 87GB 215950.2+5031417 was 
also identified a th e cata log of bl azars likely to be detected 
by Fermi f#657: [Sowards-Emmerd et alihooBh . 


5 DISCUSSION 

We have analyzed the AGILE gamma-ray point source, 
AGL J2241+4454, to search for evidence that the emission 
originates from the Be-BH binary system HD 215227. We 
find a gamma-ray flux of (1.8 ± 0.7) x 10 -6 ph cm -2 s _1 at 
3cr significance, in agreement with the initial report. Despite 
the large AGILE flux, we find no evidence of emission in the 
Fermi data using three different time binning variations. 

There are two main reasons that may be partially or 
entirely responsible for the Fermi non-detection. The first 
is that the flare was produced mostly lower energy photons 
(< a few hundred MeV) with a steep power law. The pho¬ 
tons detected by AGILE were mostly in the 100 MeV range 
where the performance of the LAT suffers from a lower effec¬ 
tive area and larger PSF compared to energies higher than 
1 GeV. Second, the flare may have been short (less than 
~6 hours) and peaked in between Fermi GTIs. This theory 
is supported by Figure [4] which shows that the peak de¬ 
tected by AGILE occurred outside of the Fermi GTIs for 
the region around the flare. Considering these possibilities, it 
seems that some combination of the first two is most likely. 
This issue may be clarified when the Fermi science team 
releases the new Pass 8 data in mid-2015. The improved 

3 http://simbad.u-strasbg.fr/simbad/ 


event reconstruction will provide significant increases in the 
effective area and field of view, especially at energies below 
100 MeV. An updated analysis with the Pass 8 data may at 
least rule out the possibility that the photons were of too 
low energy for Fermi to detect. 

In any case, it appears that the proposed gamma- 
ray source AGL J2241+4454 does not fit within the cat¬ 
egory of known gamma-ray binaries such as LS 5039 or 
LS I +61 303. The five known gamma-ray binaries all show 
orbital phase modulation of their high-ener gy (GeV) a nd 
very-high energy (TeV) gamma-ray emission ( Dubusll2013l h 
while AGL J2241+4454 does not emit gamma-rays at a con¬ 
tinuously detectable level, if at all. Instead, it is more likely 
that the HD 215227 binary system is a HMXB and micro¬ 
quasar that has not been discovered earlier owing to its 
seemingly persistent quiescent state. Deeper X-ray obser¬ 
vations will be of future importance to determine whether 
or not an ADAF model accurately describes HD 215227 or 
if other emission mechanisms are at play. Additionally, con¬ 
temporaneous X-ray and radio observations may be able to 
further test disk-jet coupling theories for the system. 

During the investigation into AGL J2241+4454, we 
identified a gamma-ray point source that was not present 
in the 2FGL catalog, but now appears in the latest 3FGL 
catalog. The source, 3FGL J2201.7+5047, is variable on 
monthly (and perhaps shorter) timescales and has been 
bright (> 10 -8 ph cm -2 s -1 ) for the past few years. We, 
and independently, the Fermi LAT team, identify the flat- 
spectrum radio quasar 87GB 215950.2+503417 as the most 
probable source of the gamma-ray emission because of its 
proximity to the gamma-ray peak, its known X-ray and radio 
emission, and its prior identification as a potential gamma- 
ray emitter. 


ACKNOWLEDGMENTS 

The Fermi LAT Collaboration acknowledges generous ongo¬ 
ing support from a number of agencies and institutes that 
have supported both the development and the operation 
of the LAT as well as scientific data analysis. These in¬ 
clude the National Aeronautics and Space Administration 
and the Department of Energy in the United States, the 
Commissariat a lEnergie Atomique and the Centre National 
de la Recherche Scientifique/Institut National de Physique 
Nucleaire et de Physique des Particules in France, the Agen- 
zia Spaziale Italiana and the Istituto Nazionale di Fisica Nu- 
cleare in Italy, the Ministry of Education, Culture, Sports, 
Science and Technology (MEXT),High Energy Accelerator 
Research Organization (KEK) and Japan Aerospace Explo¬ 
ration Agency (JAXA) in Japan, and the K. A. Wallenberg 
Foundation, the Swedish Research Council and the Swedish 
National Space Board in Sweden. We are grateful for sup¬ 
port from the National Science Foundation through grant 
AST-1109247. This research has made use of the SIMBAD 
database, operated at CDS, Strasbourg, France. Part of this 
work is based on archival data, software, or online services 
provided by the ASI SCIENCE DATA CENTER (ASDC). 
We also thank the anonymous referees for many helpful com¬ 
ments. 


© 2014 RAS, MNRAS 000, □[6] 











6 Michael J. Alexander and M. Virginia McSwain 


REFERENCES 

Abdo A. A. et al., 2009a, ApJ, 701, L123 
Abdo A. A., et al., 2009b, Science, 326, 1512 
Atwood W. B., et al., 2009, ApJ, 697, 1071 
Bodaghee A., Tomsick J. A., Pottschmidt K., Rodriguez 
J., Wilms J., Pooley G. G., 2013, ApJ, 775, 98 
Bulgarelli A. et al., 2010, The Astronomer’s Telegram, 
2512, 1 

Bulgarelli A., et al., 2011, The Astronomer’s Telegram, 
3151, 1 

Bulgarelli A., et al., 2012, A&A, 538, A63 
Casares J., Negueruela I., Ribo M., Ribas I., Paredes J. M., 
Herrero A., Simon-Di'az S., 2014, Nature, 505, 378 
Corbel S., Coriat M., Brocksopp C., Tzioumis A. K., 
Fender R. P., Tomsick J. A., Buxton M. M, Bailyn C. D., 
2013, MNRAS, 428, 2500 
Dubus G., 2013, A&ARv, 21, 64 

Gallo E., Miller B. P., Fender R., 2012, MNRAS, 423, 590 
Garcia M. R., McClintock J. E., Narayan R., Callanan P., 
Barret D., & Murray S. S., 2001, ApJ, 553, L47 
Gregory P. C., Condon J. J., 1991, ApJS, 75, 1011 
Hadasch et al., 2012, ApJ, 749, 54 

Lucarelli F. et al., 2010, The Astronomer’s Telegram, 2761, 
1 

Malyshev D., Zdziarski A. A., Chernyakova M., 2013, MN¬ 
RAS, 434, 2380 

McClintock J. E., Remillard R. A., 1986, ApJ, 308, 110 
Merrill P. W., Burwell C. G., 1943, ApJ, 98, 153 
Mori M., Kawachi A., Nagataki S., Naito T., 2013, ArXiv 
e-prints 

Munar-Adrover P., Paredes J. M., Ribo M., Iwasawa K., 
Zabalza V., Casares J., 2014, ArXiv e-prints 
Orosz J. A., McClintock J. E., Aufdenberg J. P., Remillard 
R. A., Reid M. J., Narayan R., Gou L., 2011, ApJ, 742, 
84 

Owocki S. P., Romero G. E., Townsend R. H. D., Araudo 
A. T., 2009, ApJ, 696, 690 

Remillard R. A., McClintock J. E., 2006, ARA&A, 44, 49 
Reynolds M. T., Reis R. C., Miller J. M., Cackett, E. M., 
& Degenaar, N., 2014, MNRAS, 441, 3656 
Shrader C. R., Titarchuk L., Shaposhnikov N., 2010, ApJ, 
718, 488 

Sowards-Emmerd D., Romani R. W., Michelson P. F., 
Healey S. E., Nolan P. L., 2005, ApJ, 626, 95 
Tavani M. et al., 2009, A&A, 502, 995 
van Kerkwijk M. H., Charles P. A., Geballe T. R., King 
D. L., Miley G. K., Molnar L. A., van den Heuvel E. P. J., 
van der Klis M., van Paradijs J., 1992, Nature, 355, 703 
Voges W. et al., 1999, A&A, 349, 389 
Williams S. J., Gies D. R., Matson R. A., Touhami Y., 
GrundstromE. D., Huang W., McSwain M. V., 2010, ApJ, 
723, L93 

Zdziarski A. A., Mikolajewska J., Belczyriski K., 2013, MN¬ 
RAS, 429, L104 

Zhang S.-N., 2013, Frontiers of Physics, 8, 630 


© 2014 RAS, MNRAS 000, HH] 


